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Abstract 

We demonstrate how to obtain from the Schwinger-Dyson equations of QCD an infrared finite 
ghost propagator in the Feynman gauge. The key ingredient in this construction is the longitudinal 
form factor of the non-perturbative gluon-ghost vertex, which, contrary to what happens in the 
Landau gauge, contributes non-trivially to the gap equation of the ghost. The detailed study of 
the corresponding vertex equation reveals that in the presence of a dynamical infrared cutoff this 
form factor remains finite in the limit of vanishing ghost momentum. This, in turn, allows the 
ghost self-energy to reach a finite value in the infrared, without having to assume any additional 
properties for the gluon-ghost vertex, such as the presence of massless poles. The implications of 
this result and possible future directions are briefly outlined. 

PACS numbers: 12.38.Lg, 12.38.Aw 
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I. INTRODUCTION 



The non-perturbative properties of the basic Green's functions of QCD have been the focal 
point of intensive scrutiny in recent years, with particular emphasis on the propagators of the 
fundamental degrees of freedom, gluons, quarks, and ghosts. Even though it is well-known 
that these quantities are not physical, since they depend on the gauge-fixing scheme and 
parameters used to quantize the theory, it is generally accepted that reliable information 
on their non-perturbative structure is essential for unraveling the infrared (IR) dynamics of 
QCD. 

There are two main tools usually employed in this search: the lattice, where space-time 
is discretized and the quantities of interest are evaluated numerically , and the 

intrinsically non-perturbative equations governingthe dynamics of the Green's functions, 
known as Schwinger- Dyson equations (SDE) 0,0,16, 3]- In principle, the lattice includes all 
non-perturbative features and no approximations are employed at the level of the theory. In 
practice, the main limitations appear when attempting to extrapolate the results obtained 
with finite lattice volume to the continuous space-time limit. On the other hand, the main 
difficulty with the SDE has to do with the need to devise a self-consistent truncation scheme 
that preserves crucial field-theoretic properties, such as the transversality of the gluon self- 
energy, known to be valid both perturbatively and non-perturbatively, as a consequence of 
the BRST symmetry ^. 

Significant progress has been accomplished on this last issue due to the development of 
the truncation scheme that is based on the all-order correspondence jo] between the pinch 



technique (PT) 10|, [ill] and the Feynman gauge of the Background Field Method (BFM) jl2| . 



One of its most powerful features is the special way in which the transversality of the gluon 
self-energy is realized. Specifically, by virtue of the Abelian-like Wis satisfied by the vertices 
involved, gluonic and ghost contributions are separately transverse, within each order in 
the "dressed- loop" expansion of the SDE [iS^] for the gluon propagator. This property, in 
turn, allows for a systematic truncation of the full SDE, preserving at every step the crucial 
property of gauge invariance. 

The first approximation to the SDE of the gluon propagator involves the one-loop dressed 
gluonic graphs only, since in this scheme the ghost loops may be omitted without compro- 
mising the transversality of the answer. As is well-known, the Feynman gauge of the BFM 
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is particularly privileged, being dynamically singled out as the gauge that directly encom- 



passes the relevant gauge cancellations of the PT 



9|. Therefore, the aforementioned one- loop 



dressed graphs have been considered in this particular gauge. The detailed study of the re- 
sulting integral equation for the gluon propagator gave rise to solutions that reach a finite 
value in the deep IR Q, [l^. Following Cornwall's orig inal idea jiol. [isl of describing the IR 



sector of QCD in terms of an effective gluon mass [16|, [ITJ, these solutions have been fitted 
using "massive" propagators of the form A~^(g^) = + m^(g^), with m^(0) > 0, and the 
crucial characteristic that m?{q^) is not "hard", but depends non-trivially on the momentum 
transfer q^. In addition, finite solutions for the gluon propagator in the Landau gauge have 
been reported in various lattice studies [l8| , and were recently confirmed using lattices with 
significantly larger volumes fiol . 

Even though the omission of the ghost loops within this formulation does not introduce 
any artifacts, such as the loss of transversality, the actual behavior of the ghosts may change 
the initial prediction for the gluon propagator, not just quantitatively but also qualitatively. 
For example, an IR divergent solution for the ghost propagator could destabilize the finite 
solutions found for the gluon propagator. Therefore, a detailed study of the ghost sector 
constitutes the next challenge in this approach. In the present work we will consider the SDK 
for the ghost sector in the (BFM) Feynman gauge, in order to complement the corresponding 
analysis presented in 1^, 1^ in the same gauge. The BFM Feynman rules are in general 
different to those of the covariant renormalizable gauges 12|; in the former, for example, 
in addition to the bare gluon propagator, the bare three- and four-gluon vertices involving 
background and quantum gluons depend on the (quantum) gauge fixing parameter. Notice, 
however, that, since there are no background ghosts, the Feynman rules relevant for the ghost 
sector are identical to both the covariant gauges and the BFM. Therefore, the analysis and 
the results presented in this article carries over directly to the conventional Feynman gauge. 

In this article we demonstrate that the ghost propagator in the Feynman gauge can be 
made finite in the IR, through the self-consistent treatment of the gluon-ghost vertex and 
the ghost gap equations. The key ingredient that makes this possible is the "longitudinal" 
form-factor in the tensorial decomposition of the gluon-ghost vertex, W^^'^ij), q, k), i.e. the co- 
factor of kfj_, where k is the four-momentum of the gluon; evidently this term gets annihilated 
when contracted with the usual transverse projection operator. As we will explain in detail, 
this component acquires a special role for all values of the gauge fixing parameter, with the 
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very characteristic exception of the Landau gauge. The reason is simply that in the Landau 
gauge the entire gluon propagator is transverse, both its self-energy and its free part, whereas 
for any other value of the gauge-fixing parameter the free part is not transverse. As a result, 
when the gluon-ghost vertex is inserted into the SDE for the ghost propagator, D{p'^), its 
part proportional to dies when contracted with the gluon propagator in the Landau gauge; 
however, in any other gauge it survives due to the free-part of the gluon propagator. The 
resulting contribution has the additional crucial property of not vanishing as the external 
momentum of the ghost goes to zero. Therefore, contrary to what happens in the Landau 
gauge where only the part of the vertex proportional to survives, one does not need to 
assume the presence of massless pole terms of the form in order to obtain a nonvanishing 
value for D~^(0). Instead, the only requirement is that the longitudinal form factor simply 
does not vanish in that limit. 

The paper is organized as follows: In section [TTl we set up the SDE for the ghost propa- 
gator, assuming the most general Lorentz structure for the fully dressed gluon-ghost vertex 
r|^^'^(p, g, fc). We then discuss under what condition the resulting expression may yield a 
finite value for D~^(0), and analyze the profound differences between the Landau and the 
Feynman-type of gauges. In section IIIII we first derive the gluon-ghost vertex under certain 
simplifying assumptions, discuss in detail the approximations employed. Next we study its 
non-perturbative solutions employing various physically motivated, IR finite Ansatze IR for 
the gluon and ghost propagators. In Section [IV] we combine the results of the previous 
two sections, deriving the self-consistency condition necessary for the system of equations 
to be simultaneously satisfied. Finally, in section |V] we discuss our results and present our 
conclusions. 

II. GENERAL CONSIDERATIONS ON THE IR BEHAVIOR OF THE GHOST 

In this section we derive the SDE for the ghost propagator D{p^) in a general covariant 
gauge, and study qualitatively its predictions for -D(O) for various gauge choices. In particu- 
lar, we establish that away from the Landau gauge the ghost propagator may acquire a finite 
value at the origin, without the need to assume a singular IR behavior for the form factors 
of the fully dressed ghost-gluon vertex entering into the SDE. Our attention will eventually 
focus on the Feynman gauge, which, as mentioned in the Introduction, is singled out within 
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the PT-BFM scheme. 

The full ghost propagator D"-^[p) is usually written in the form 

D''\p) = iS^^Dip) , (2.1) 

and the SDE satisfied by DijP'), depicted diagrammatically in FigUl reads 

D-^p") =p' + iCaQ^ j [dk] V^^^^,{k)T''{p,p + k, k)D{p + k) . (2.2) 

We have used f°-'^'^f'"^'^ = S°'^Ca, with Ca the Casimir eigenvalue in the adjoint representation 
[Ca = N for SU (N)], and have introduced the short-hand notation [dk] = d'^k/ {2nY , where 
d = 4 — e is the dimension of space-time used in dimensional regularization. Af^^{k) is the 
fully dressed gluon propagator, whereas F denotes the fully dressed gluon-ghost vertex and 
r its tree-level value. 



k 




FIG. 1: The SDE of the ghost propagator. 



Specifically, in the covariant gauges the full gluon propagator A^{,(A;) = —i5'^^A^y{k) has 
the general form 



P,.{k)A{k^)+^^ 



where 



k ky 

P^iu{k) = g^y 



k 



2 ' 



(2.3) 



(2.4) 



is the transverse projector, and ^ is the gauge fixing parameter; = 1 corresponds to the 
Feynman gauge and = to the Landau gauge. The scalar function A(A;^) is related to the 
all-order gluon self-energy Yi^y{k), 



through 



n^.(A;) = V^,{k)Tl{e) , 



(2.5) 



(2.6) 



The bare gluon-ghost vertex appearing in (12. 2p is given by T^^-^ = —gf^"'^qfj,, with 
{q = p + k). Choosing and A;^ as the two hnearly independent four- vectors, the most gen- 
eral decomposition for the fully dressed gluon-ghost vertex T^'^'^ij), q, k) is expressed as [20| 



K^P^q.k) = -gr^'T.ip^q.k), 
T^{p,q,k) = A{p\q',e)p^ + B{p',q',e)k 



(2.7) 



where k is the outgoing gluon momentum, and p, q the outgoing and incoming ghost mo- 
menta, respectively. The dimensionless scalar functions A{p'^ , q'^ , k'^) and B{p'^,q'^,k'^) are 
the form factors of the gluon-ghost vertex. In particular, notice that the tree-level result 
is recovered when we set Alp"^, g^, k"^) = 1 and B{p'^, g^, k'^) = 0. Finally, it is important to 
emphasize that all fully-dressed scalar quantities {D, A, A, and B) depend explicitly (and 
non-trivially) on the value of the gauge-fixing parameter ^ already at the level of one-loop 
perturbation theory. 

It is then straightforward to derive the Euclidean version of Eq. (l2.2p : to that end, we 
set p2 = — p|, define Ae(p|) = — A(— p|), and De{Pe) = —D{—p%) , and for the integration 
measure we have [dk] = i[dk]E = id'^kE/ {^t^Y- Suppressing the subscript "E" everywhere 
except in the integration measure, and without any assumptions on the functional form of 
A(p2, A;2) and B{p^, g^, k"^)^ the ghost SDE of Eq. becomes 



D-\p') = p'-C^g'Jldk] 
-CAg'd [dk] 



p 



{p ■ kf 

fc2 



A{p\q\k')A{k)D{p + k) 



p ■ k 



A{p\q\e) + B{p\q\k') + ^-^ 



fc2 



D{p + k) 



-CAg^ij[dk]^B{p\q\e)D{p + k) 



(2.8) 



As a check, we can recover from (12. 8p the one-loop result for the ghost propagator in the 
Feynman gauge (^ = 1) by substituting the tree- level expressions for the ghost and gluon 
propagators and setting AijP', q^, k"^) = 1 and B{p'^, q"^, k"^) = 0; specifically 



D-\p'' 



p 



1 + ^lnf^ 



327r2 



(2.9) 



In order to obtain from (12. Sp the behavior of D{jP') for the full range of the momen- 
tum p2 one needs to provide additional information for the forms factors A{p'^^ q^, k"^) and 
B{p'^,q'^,k'^), obtained from the corresponding SDE satisfied by the gluon-ghost vertex. 



Thus, the complete treatment of this problem would require the solution of a complicated 
system of coupled SDE. However, several interesting conclusions about the IR behavior of 
D{p^) may be drawn, by considering the qualitative behavior of the forms factors A{p'^, g^, k"^) 
and g^, fc^), as p ^ 0. 

We start by considering what happens in the Landau gauge. First of all, let us assume 
that the various quantities appearing on the r.h.s. of (12.81) are regular functions of ^ 21|. 
Then, if we set ^ = 0, only the first integral on the r.h.s. of (12.81) survives; thus, D~^{p'^) is 
only affected by the functional form of g^, k"^). In particular, the behavior of D{p^) as 
p will depend on whether g^, fc^) is divergent or finite in that limit, i.e. on whether 
or not A{p'^,q'^,k^) contains terms. Evidently, if A{p'^,q'^,k'^) does not contain poles, 

one has that lim D^^{0) = 0, and therefore the ghost propagator will be divergent in the IR. 
On the other hand, if Alp"^ , q"^ , k"^) contains terms, limD~^(0) 7^ allowing for finite 

solutions for the ghost propagator. 

According to this general argument, the only way for getting an IR-finite propagator 
in the Landau gauge is by assuming that y4(p^,g^,fc^) contains poles (22|, l23|. However, 
lattice simulations in the Landau gauge seem to favor a IR-finite A{p'^, g^, k'^); specifically, it 
was found that deviations of the gluo n-ghost vertex from its tree-level value are very small 



in the IR, i.e. g^, fc^) ^ 1 [2J]. In addition, a detailed study of the SDE equation 

for r in the same gauge shows no singular behavior for y4(p^,g^,fc^) 25|] . These findings 
appear to be consistent with recent lattice results on the non-perturbative structure of the 
ghost propagator, which indicate that D^^ijP') in the Landau gauge diverges, at a rate that 
deviates only mildly from the tree-level expectation of l/p^ 19[ |. 

Evidently, the picture for ^ 7^ is drastically different. Indeed, away from the Landau 
gauge the r.h.s of (12. 8p involves both form factors, g^, fc^) and -B(p^, g^, k"^). Moreover, 
unlike the first two terms, the third one does not contain any kinematic factors proportional 
to p. Thus, in order for it not to vanish as p ^ one does not need to assume any singular 
structure for _B(p^, g^, /c^); instead, it is sufficient to simply have that 5(0, fc^, fc^) 7^ 0. 

After this key observation, we will take the limit of of Eq. (12.81) as p — > 0, assuming 
that A(p^, g^, fc^) does not contain (1/p^) terms. Focusing for concreteness on the physically 
relevant case of ^ = 1, we find that in the aforementioned kinematic limit Eq. (l2.8p reduces 
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to 



D-\Q) = -Ca9 J [dk]^ 5(0, r , D{k) . (2.10) 

Of course, if the assumption that A{p'^ , q"^ , k'^) is regular as p — > does not hold, then the 
other integrals will also contribute to the r.h.s. of fl2.10p . However, modulo the rather 
contrived scenario of fine-tuned cancellations, the r.h.s. will still be different from zero. 
Evidently, from fl2J0D we deduce that if 5(0, fc^, A;^) = than D-^{0) = 0. On the other 
hand, if 5(0, ^ 0, i.e. if it does not vanish identically, then one may have a non- 

vanishing D^^{0). Of course, having a non-vanishing 5(0, k"^, k^) is not a sufficient condition 
for 5^^(0) 7^ 0; one has to assume in addition that (i) the integral on the r.h.s. of (12.101) . 
is convergent, or it can be made convergent through proper regularization, and (ii) that 
the integral is not zero due to some other, rather contrived circumstances (for instance, 
if 5(0, A;^, A;^) turned out not to be a monotonic function, the various contributions from 
different integration regions could cancel against each other). 



An explicit calculation may confirm that 5(0, k"^, k"^) vanishes at one-loop 26|, and it is 
reasonable to expect this to persist to all orders in perturbation theory. Therefore, in what 
follows we will examine the possibility that 5(0, k"^, k"^) may not vanish non-perturbatively. 
In particular, we will study the SDE determining B{p'^,q'^,k'^) for the special kinematic 
configuration appearing in fl2.10p . namely where the outgoing ghost momentum, p, is set 
equal to zero (i.e. p = and q = k). In the context of the linearized approximation that 
we employ in the next section this kinematic configuration offers the particular technical 
advantage of dealing with a function of only one variable instead of two. 

III. THE GLUON-GHOST VERTEX 

In this section we set up and solve, after certain simplifying approximations, the SDE 
governing the behavior of the form factor 5(0, /c^, A;^). This can be done by taking the 
following limit of the gluon-ghost vertex, r^(p, q, k) , 

" 1 



B{0,k\k') = lim 



^^kf'T^ip^q^k) 



(3.1) 



where r^(p, q, k) obeys the SDE Q represented in FigE] 

We next introduce some approximations regarding the form of the two-ghost-two-gluon 
scattering kernel, appearing on the r.h.s. of Fig. ([2]). The first approximation is to keep only 




FIG. 2: SDE for the gluon-ghost vertex. 

the lowest order contributions in its skeleton expansion, i.e. we expand the aforementioned 
kernel in terms of the IPI fully dressed three-particle vertices of the theory, neglecting 
diagrams that contain four-point functions. 

We then arrive at the truncated SDE shown in Fig. ([3]), which reads. 




(fli) (02) 



FIG. 3: Truncated version of the SDE for the gluon-ghost vertex. 

K^ip, q, k) = r';" + T'^'ip, q, k) + Tl^ip, q, k) , (3.2) 
where the closed expressions corresponding to the diagrams (oi) and (02) are given by 
K"" la. = / \dl] T^:r\l + P,l + q, k)D,e'{l + P) Ki-\p, l+P,l) KnAl) C'^" + q) , 

r^"^ L, = j [di] rfr q-i,i- p)AzUi - p) ir^?'"^'(p, i.i-p) Dnwii) ry'A:::{i - q) , 

(3.3) 

with the momentum routing as given in Fig. ([3]). 

Our next approximation is to linearize the equation by substituting in (13. 3p 
rj;'"'^(/ + p^l + q^k) and W^J;^{—k, q — l,l—p) by their bare, tree-level expressions. Since we 
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are eventually interested in the limit of the equation as j9 ^ 0, this amounts finally to the 
replacement 



K^i-k, q-l.l-p) ^ gr'^m - k),g,^ - (A; + lU,, + {2k - lU,,] . (3.4) 

in diagrams (ai) and (02), respectively. The diagrammatic representation of the resulting 
contributions at p — * is given in Fig. (jl]) . 



>T k 



l + k 



/ \ 

/ 



\ 



/ 

/ 

/ c 



*x/\/\/\0/\/\/\/^ 




(«2) 



FIG. 4: Contributions for the gluon-ghost vertex equation in the hmit of p — > 0. 



Factoring out the color structure by using the standard identity f°-^"^fbm'>T'f'^"'^ = 1 CaJ"''"^, 



it is easy to verify that in the limit p the linearized version of Eg. (13. 3 p reads 



irr(o, k, k) 



a2 



bed 



Ca9' 



[dl] + fc),,/,A""'(/)S(0, /^ l^)D(l)D{l + k) , 
j [dl] T^,, U,>^'"'' (/) A^^' (/ - A;)5(0, l\ l^)D{l) . (3.5) 



Since the bare gluon-ghost is proportional to p^, it follows immediately from Eqs. (l3.1l) . 
and ([33D, that 



k^' 



B{0,k\k') = — T^iO,k,k) l^+T^{0,k,k) 



W^(0,fc,fc) = --Ca9' [dl] 



W^(0,fc,fc) 



k-l + 
\k ■ If 



{k-lf 



B{0,f,f)D{l)D{l + k), 
B{0,f,P)D{l)A{l + k). (3.6) 
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The Euclidean version of (13.60 can be easily derived using the same rules as before, leading 
to 

+ jd^lB{0,f,l^)D{l)A{l + k) + ^ JdH {k ■ l)B{0,l^,l^)D{l)D{l + k)"^ .(3.7) 

It is convenient to express the measure in spherical coordinates, 

dH = 27T dxsin'x / dyy; (3.8) 
Jo Jo 

and rewrite (13.71) in terms of the new variables x = k'^, y = P, and z = {I + k)'^. In order to 
convert Eq. (13.71) into a one- dimensional integral equation, we resort to the standard angular 
approximation, defined as 



eix-y)fix) + 9iy-x)fiy) 



(3.9) 



where d{x) is the Heaviside step function. 

Then, introducing the above change of variables and using Eq. (13.81) and (13. 9p in (13.71) . we 
arrive at the following linear and homogeneous equation 

B{0,x,x) = ^^\^[D{x)-A{x)] rdyy'B{0,y,y)D{y) 
1287r^ [x^ -'Jo 

X) /"OO 

dy{x-2y)BiO,y,y)Diy)[Diy)-Aiy)]+2 / dyyBiO,y,y)Diy)Aiy) 



--D{x) I dyy^B{0,y,y)D{y) + 4A{x) I dyyB{0,y,y)D{y) } . (3.10) 



X 

Due to the linear nature of (13.101) it is evident that if B is one solution then the entire 
family of functions cB, generated by multiplying B by an arbitrary constant c, are also 
solutions. 

Before embarking into the numerical treatment of (]3.10p . it is useful to study the asymp- 
totic solution that this equation furnishes for x — oo. In this limit on can safely replace the 
various propagators appearing on the r.h.s of (I3.10p by their tree-level values, i.e. A{t) 1/t 
and D(t) 1/t with (t = x,y). Then, the first and second terms vanish, and the leading 
contribution comes from the third term of (I3.10p . Specifically, the asymptotic behavior of 
B{0,x,x) is determined from the integral equation 

B{0,x,x) = X Tdy^^^^^, (3.11) 

Jx y 
11 



where A = C ^g"^ / QAti"^ . Eg. (13. Ill) can be solved easily by converting it into a first-order 
differential equation, which leads to the following asymptotic behavior 



-8(0, X, x) = ax 



(3.12) 



with a is an arbitrary parameter, with dimension [M^]^, where M is an arbitrary mass-scale. 
As we will see in what follows, a will be treated as an adjustable parameter, whose dimen- 
sionality will be eventually saturated by that of the effective gluon mass, or, equivalently, 
by the QCD mass scale A. 

With the asymptotic behavior (13.121) at hand, we can solve numerically the integral 
equation given in (I3.10p . To do so, we start by specifying the expressions we will use for the 
gluon and ghost propagators. 

As has been advocated in a series of studies based on a variety of approaches, the gluon 

This type of behavior has been 



propagator reaches a finite value in the deep IR [27 



observed in Landau gauge in previous lattice studies 



and more recently in new, large- 
191 . Within the gauge- invariant truncation scheme implemented by 
the PT, the gluon propagator (effectively in the background Feynman gauge) was shown to 



volume simulations 



saturate in the deep IR 
a propagator of the form 



13 



14l |. The numerical solutions may be fitted very accurately by 
1 



(3.13) 



^2 _|_ ^2^^2-) ' 

where m?{k'^) acts as an effective gluon mass, presenting a non-trivial dependence on the 
momentum k'^. Specifically, the mass displays either a logarithmic running 

-1-71 



m\k' 



■rrir, 



In 



k"^ + pml 
A2 



In 



A2 



(3.14) 



where 71 > is the anomalous dimension of the effective mass, or power-law running of the 
form 



m\k' 



rrin 



rrin 



In 



k"^ + pml 



In 



72-1 



(3.15) 



A2 J / V A2 

with 72 > 1. Which of these two behaviors will be realized is a delicate dynamical prob- 
lem, and depends, among other things, on the specific form of the full three-gluon vertex 
employed in the SDE for the gluon propagator (for a detailed discussion see [ij])- Here we 
will employ both functional forms, and study the numerical impact they may have on the 
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solutions of (13.101) . A plethora of plienomenological studies favor values of mo in the range 
of 0.5-0.7GeV. 

In addition, when solving (13.101) an appropriate Ansatz for the ghost propagator D{k^) 
must also be furnished, given that we are in no position to solve the ghost SDE of (12. 8p for 
arbitrary values of the momentum, since this would require the solution of a coupled system 
of several integral equations involving D, A, and B, for arbitrary values of the four-momenta. 
Given that our aim is to study the self-consistent realization of an IR finite ghost-propagator, 
it is natural to employ an Ansatz in close analogy to (I3.13p . namely 

1 



(3.16) 



where M'^{k'^) stands for a dynamically generated, effective "ghost mass". Evidently, 
L)-i(O) = M2(0), and D-^{0) ^ provided that M^{0) ^ 0. Of course, once the corre- 
sponding solutions for B{0,x,x) have been obtained the self-consistency of the Ansatz for 
M'^{k'^) must be verified. The way this will be done in the next section is by substitut- 
ing B{0,x,x) into the (properly regularized) integral on the r.h.s. of Eq. (I2.10p . and then 
demanding that its value is equal to the M'^{0) appearing on the l.h.s. 

For the actual momentum dependence of the effective ghost mass, M{k'^) we will assume 
three different characteristic behaviors and will analyze the sensitivity of -8(0, x, x) on them. 

We will employ the following three types of M{k^): 

(i) "hard mass", i.e. a constant mass with no running, 



; 



(3.17) 



(ii) logarithmic running of the form 



M\k') = M( 



In 



A2 



A2 



— 1 — Kl 



(3.18) 



(iii) power-law running, given by 



Mo' 



k^ + Ml 



In 



e + pMi\ /^^fpMi 



A2 



A2 



K2 — 1 



(3.19) 



Clearly, the last two possibilities, (13.181) and (13.191) . are exactly analogous to the correspond- 
ing two types of running of the gluon mass, (13.140 and (I3.15p . respectively. 
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We then solve numerically Eg. (13. 101) using the gluon and ghost propagators given by 
Eqs.f l3.13|) and (13.161) . respectively, supplemented by the various types of running for m^(fc^) 
and M^(A;^). The integration range is split in two regions, [0, s] and (s, oo], where s ^ A^. 
For the second interval we impose the asymptotic behavior of (I3.12p . choosing a value for a. 

It turns out that the numerical solution obtained for B{0 •j X ^ X ^ IS rather insensitive to the 
form of the gluon mass employed, and it mainly depends on the form of the ghost propagator. 
More specifically, we can fit the numerical solution with an impressive accuracy by means 
of the simple, physically motivated function 

^(0, X, x) = ^prr^ , (3.20) 

regardless of the form of momentum dependence employed for M^(x). Evidently, for large 
values of x the above expression goes over the asymptotic solution of Eq. (l3.12l) . In Figs.(l5l), 
we present a typical solution for i?(0,a;,a;) together with the fit given by (13.201) . 




FIG. 5: The black solid line is the numerical solution of Eq. (j3.10p . assuming logarithmic type of 
running for m?{k'^) and M^{P), with 71 = ki = 0.6, = 0.35 GeV^, = 0.4 GeV^, p = 4, and 
o"^/^ = IGeV^. The red dashed line represents the fit of Eq. (j3.20p : the relative difference between 
the two curves is less than 1% (note the fine spacing of the y axis). 
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IV. INFRARED FINITE GHOST PROPAGATOR 



In the previous section we have obtained the general solutions for B{0,x,x), under the 
assumption that the ghost propagator was finite in the IR, and more specifically that it 
was given by the general form of fl3.16p . The next crucial step consists in substituting the 
solutions obtained for B{0,x,x) into fl2.10p and in examining under what conditions the 
two hand sides of the equation can be made to be equal. As we will see this procedure will 
eventually boil down to constraints on the values that one is allowed to choose for the free 
parameter a. 

Substituting Eqs. fl3.16l) and fl3.20l) into fl2.1Up . we arrive at 

1 



D-\0) = -Ca9'cj J [dk] 
The r.h.s. of (14.21 1) is simply given by 



2^ll+A 



(4.21) 



D-\0) = M( 



5 



(4.22) 



for any form of M^(A;^). Let us first verify the self-consistency of 04.211) for the case where 
the ghost mass vanishes identically, i.e. M^(fc^) = 0. Then, (I4.2ip reduces to nothing but 
the standard dimensional regularization result 29| 



[dk]{k 



2\-a 



0. 



(4.23) 



valid for any value of a, for the special value a = 1 + A. 

For non- vanishing M'^{k'^) the integral on the r.h.s. of (14.2 ID is UV divergent: at large k'^ 
it goes as (Auv)^''"'^, where Auv is a UV momentum cutoff. It turns out that the r.h.s. can 
be made UV finite by simply subtracting from it its perturbative value, i.e. the vanishing 
integral of (1123]) [30 1. 

Carrying out this regularization procedure explicitly, one obtains 

1 1 



= ~CA9^a J [dk] 



[P + M2(P)]i+^ 

[dk] 



[A;2 + M2(P)] 



M2(A;2 



k^ 



(4.24) 



It is now elementary to verify that the integral on the r.h.s of ( 14.24^ converges. At large 
k"^ we can expand the second term in the parenthesis and neglecting in the denominator 
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M'^{k'^) next to k'^, we find that the resulting integral (apart of multiplicative factors) is 
given by 

dy-^. (4.25) 

Notice that the above integral converges even for the less favorable case of a constant M^(?/); 
then, f l4.25p is proportional to ?/~^, and is therefore convergent, since A > 0. Clearly, when 
M^(?/) drops off in the UV, as described by f l3.18p or (13.191) . the integral converges even 
faster. Next we will analyze separately what happens for each one of the three different 
Ansatze we have employed for M^(?/), Eqs. (l3.17l) - (I3.19p . 

The case of a constant ghost mass can be easily worked out. Replacing M^(A;^) Mq 
in Eg. (14.241) . keeping only the leading contribution to the integral, we arrive at (notice the 
cancellation of the coupling constant appearing in front of the integral) 

= J^Mo'(^-") . (4.26) 
1 — A 

Then, in order to enforce the equality of both sides of (14.261) a must satisfy 

a = il—^M^^ . (4.27) 

Evidently, a depends very weakly on Mq, and its value is practically fixed at 1/4. Indeed, 
given that A is a small number, of the order of O(10~^), Eq. (l4.27p may be expanded as 



(4.28) 



from where it is clear that a can only assume values slightly different of 1/4. In Fig. ([6]), we 
show this mild dependence of a on Mq, for A = 300 MeV. 

We next turn to the case where M'^{y) displays the logarithmic or power-law dependence 
on the momentum, described by Eqs. (I3.18P and (I3.19p . respectively. Now the integrals 
cannot be carried out analytically and have been computed numerically. Choosing different 
values for ki, K2, and p, we obtain the curves presented in Fig. ([7]) and Fig. ([8]), showing the 
dependence of a on Mq. 
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FIG. 6: a as a function of the hard ghost mass Mq, obtained from Eq. (|4.27p . 



Several observations are in order: 

(i) For both types of running the results show a stronger dependence on Mq than in the 
case of the hard mass. 

(ii) The range of possible values for a increases significantly. Whereas in the case of 
constant mass one was practically restricted to a unique value for a, namely cr ^ 1/4 [viz. 
Fig-®], now one may obtain self-consistent solutions choosing values for cr over a much 
wider interval. 




0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 



M„ [GeV] 

FIG. 7: a as function of Mq, when M'^{k'^) runs logarithmically, as in Eq. ()3.18p . 
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0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 



M„ [GeV] 

FIG. 8: a as function of Mq, when the power-law running of Eg. ([31^ is assumed for M'^{k^). 

(ill) There is a qualitative difference between the logarithmic and power-law running: 
in the former case a is a decreasing function of Mq, while in the latter it is increasing. 
This offers the particularly interesting possibility of finding values for a that furnish self- 
consistent solutions for either types of running of M^{k^). A characteristic example where 
Eg. (14.241) is satisfied for the same value of Mq for both types of running is shown in Fig. ({9]): 
for a ^ 20 one may generate a ghost mass of Mq ^ 560 MeV, assuming for M'^{k'^) either 
the logarithmic running of Eq. fl3.18p . or the power-law running of Eq. fl3.19p . 



60- 



50 - 



- Power-law k^=2.8 and p=2 

- Logarithmic k^=1 .3 and p=2 



, 30 



20- 



0-1 , 1 , 1 , 1 , 1 , 1 , 1 , 

0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 

M„ [GeV] 



FIG. 9: For a = 20, a ghost mass can be generated from Eq. (j4.24|) for either type of running. 
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V. DISCUSSION AND CONCLUSIONS 



In this article we have demonstrated that it is possible to obtain from the SDEs of QCD 
an IR finite ghost propagator in the Feynman gauge. In this construction the longitudinal 
component of the gluon-ghost vertex, which is inert in the Landau gauge, assumes a central 
role, allowing for -D(O) to be finite. This is accomplished without having to assume any 
special properties of the form-factor, other than a nonvanishing limit in the IR; in particular, 
we do not need to impose the presence of massless poles of the type 

Our procedure may be summarized as follows. First of all, since we are interested in the 
possibility of obtaining D^^{0) 7^ we have focused on the form of the ghost gap equation 
in the limit of vanishing external momentum, p ^ 0. Next we have linearized the SDK 
for the form-factor B{p'^,q'^,k'^) , and have looked for solutions for the special kinematic 
configuration of vanishing ghost momentum, -6(0, k^, A;^), which is relevant for the ghost gap 
equation. The solution may be fitted in the entire range of momenta with a particularly 
simple, physically motivated expression. Coupling the two equations together, we have 
obtained the conditions necessary for self-consistency. It essentially boils down to relations 
between the free parameter a and the values of D^^{0), or equivalently Mg, as captured 
in Figs (I6])-(IH]). These figures furnish the value of Mq one obtains if a concrete value of a 
is chosen, assuming certain characteristic types of running for the the ghost mass function 
M^(/c^). The freedom in choosing the value of a will be restricted, or completely eliminated, 
in the non-linear version of the vertex equation. It would certainly be interesting to venture 
into such a study, because it is liable to pin down completely the value of D^^{0). 

The most immediate physical implication of the results presented here is that the finite 
gluon propagator obtained in the previous SDK studies in the PT-BFM framework, with 
the ghost contributions gauge-invariantly omitted, will not get destabilized by the inclusion 
of the ghost loops. Specifically, one would expect that the addition of the ghost loop into 
the corresponding SDK should not change the qualitative picture. The quantitative changes 
induced should also be small; mainly the correct coefficient of llC^/487r^ multiplying the 
renormalization group logarithms will be restored (without the ghosts it is 10Ca/4:8tt'^), and 
it might inflate or deflate slightly the corresponding solutions for the gluon propagator in 
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the intermediate region between 0.1-lGeV^ Of course, a complete analysis of the coupled 
SDE system is needed in order to fully corroborate this general picture. 

Given the complexity and importance of the problem at hand it would certainly be 
essential to confront these SDE results with lattice simulations of the ghost propagator in 
the Feynman gauge. In addition, since the formulation of the BFM on the lattice has been 



presented long ago by Dashen and Gross (in the Feynman gauge) |3l|], and has already been 



used 
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33( 1 . one might also consider the possibility of simulating the gluon propagator 
within that particular gauge-fixing scheme, thus enabling a direct comparison with the SDE 
results predicting an IR finite answer. 
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